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Abstract 



It is shown that field-theory based single boson exchange potentials cannot be 

identified to those of the Yukawa or Coulomb type that are currently inserted in the 

^ ' Schrodinger equation. The potential which is obtained rather corresponds to this 

1^ , current single boson exchange potential corrected for the probability that the system 

under consideration is in a two-body component, therefore missing contributions 

due to the interaction of these two bodies while bosons are exchanged. The role of 

^ . these contributions, which involve at least two boson exchanges, is examined. The 

H I conditions that allow one to recover the usual single boson exchange potential are 

given. It is shown that the present results have some relation: i) to the failure of 

the Bethe-Salpeter equation in reproducing the Dirac or Klein-Gordon equations 

in the limit where one of the constituent has a large mass, ii) to the absence of 

corrections of relative order a log — to a full calculation of the binding energy in 

the case of neutral massless bosons or iii) to large corrections of wave-functions 

calculated perturbatively in some light-front approaches. 

PACS numbers: 13.40.Fn, 13.75.Cs, 25.45.De 
Keywords: potential, boson exchange, field theory 



1 Introduction 

It is well known that the Bethe-Salpeter equation |l|, in both the ladder approximation 
and in the limit where one of the two interacting particles has a large mass, does not 
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reproduce the Dirac or Klein-Gordon equations ||^. These are important benchmarks, 
as they accurately describe the interaction of charged particles with spin 1/2 (electrons, 
muons,...) or spin (pions,...) in the Coulomb field produced by a nucleus. The remedy 
to the above failure resides in the contribution of second order crossed diagrams (see refs 
0, ^ for some review). At first sight, taking into account these diagrams is somewhat 
surprising as both the Bethe-Salpeter and the Dirac or Klein-Gordon equations imply 
the iteration of what looks like a single boson exchange, allowing one to sum up the 
contribution of ladder diagrams. There is no obvious indication that these equations 
explicitly involve crossed diagrams. One can thus expect that the role of crossed diagrams 
in restoring the equivalence of the Bethe-Salpeter and Dirac or Klein-Gordon equations 
is partly fortuitous and results from some cancellation of contributions, which may not 
have a general character. 

It is also well known that the Bethe-Salpeter equation or the light-front approach, 
when applied to the Wick-Cutkosky model (two scalar particles exchanging a massless 
spin particle [^, Q , leads to corrections of the binding energy with an unusual character 
0, ^ 1^. In the non-relativistic limit, this model is equivalent to the Coulomb problem 
and the binding energy can be expressed in terms of the quantity, ma^, where a can be 
identified to the similar quantity used in QED [a is equal to the ratio of the coupling 
constant (?^, generally referred to in the model, with the quantity 4m^.47r). With these 
definitions, the corrections are of order m o? log - (relative order a log - with respect to 
the binding energy). These ones have been determined theoretically and are present 
in some numerical calculations [0, ||, 0. It has been suggested that the difference with 
the Coulomb interaction case, where these corrections are known to be absent, would be 
due to the spin 1 nature of the photon [^]. 

A third point we would like to mention deals with the determination of a relativistic 
wave function in the light-front approach for a realistic physical system. This represents 
a tremendous task, especially in the case of the NN system where the exchange of sev- 
eral mesons with different spins or isospins has to be considered. Present perturbative 
calculations using current deuteron wave functions as a zeroth order require a large renor- 



malization of the first order wave function, by about 20% [^, |T0|- This is much larger than 
the corrections due to incorporating relativistic kinematical effects, such as those due to 
factors (^)^''^, which amount to a few percent [0. 

This is also much larger than corrections due to renormalization effects, which are 
related to the meson in flight content of the deuteron state and have been found to be 
of the order of 3 — 4% W2. O]. Furthermore, the use of this wave function calculated 



perturbatively has revealed bad features such as an unacceptable ratio of the deuteron 
asymptotic normalizations A£, and As-, ^, or a deuteron quadrupole moment, both 
exceeding experiment by 20 — 25%. These discrepancies are too large to be considered 
as conventional relativistic corrections of the order (-)^. They rather suggest that the 
comparison is biased by the absence of essential ingredients in the calculation, which 
should be flrst completed to provide physically relevant results. 

Obviously, as mentioned in |T^, one expects these difficulties to be removed when the 
problem is treated in its totality, implying the determination of an appropriate nucleon- 



nucleon interaction model, as made in the full (energy dependent) Bonn model ||T^. We 



nevertheless believe that it is important to understand the origin of the large renorma- 
lization as well as of unexpected results in relation with the deuteron D-state. On the 
one hand, the difficulty in determining a nucleon-nucleon interaction model invites to rely 
as much as possible on what has already been achieved in the past. On the other, this 
understanding should avoid to blindly absorb in the parameters of a fitted single meson 
exchange potential contributions that are too specific to be treated in a so crude way. 
This remark is relevant for the 27r-exchange part. Although it concerns a case which is 
less realistic than the NN one, a relatively large correction of the wave function was also 
found in the light-front approach applied to the Wick-Cutkosky model in the limit of a 
small coupling P, ^ . This correction is located in the low momentum range of the wave 



function and, contrary to that one anticipated in [|T5[, which rather occurs in the high 
momentum domain, it does not go to with the momentum. For some part, it can be 
accounted for by a renormalization of the coupling constant, «—*>«(! — -a log -). 

All the problems briefiy sketched above have been raised in the framework of relativistic 
descriptions of a two-body system, such as the Bethe-Salpeter equation, or the light-front 



equation as used in [0, |T4|, |T5|. What we would like to show here is that they are rather 
due to the field-theory treatment which underlies them and that they have a much more 
general character, with a close relation to dynamical effects in the field of the many body 
problem and the related derivation of effective degrees of freedom. Essentially, the effects 
under consideration rather have a static character and are different from usual relativistic 
effects that are velocity dependent and represent a higher order effect. To emphasize our 
point, we will essentially show how the same problems occur and are solved in an approach 
which is mainly a non-relativistic one. While doing so, we are aware that the distinction is 
not straightforward. A field-theory treatment, which implies the creation and destruction 
of particles, is one of the ingredients that are required to fulfill the Lorentz covariance 
of some approaches like that followed in [|r^, |n|. However, as evidenced by the very 



existence of relativistic quantum mechanics approaches [0 , this property can be ensured 
without relying on field theory. On the other hand, a consequence of the present work 
should be to provide a better zeroth order approximation for perturbative calculations 
using energy dependent interaction models and aiming to study relativistic effects. 

The plan of the paper is as follows. In the second section, starting from a field-theory 
based single boson exchange contribution to the two-body interaction, we remind how one 
can derive an energy independent effective interaction to be employed in a Schrodinger 
equation. The difference with the current single boson exchange potential of the Yukawa 
or Coulomb type is emphasized. The third section is devoted to the corrections that the 
above difference produces, namely some renormalization of the two-body wave function 
and a^ log - terms in the binding energy in the case of an exchange of a massless (scalar) 
boson. In the fourth part, we consider the contribution due to including two boson 
exchanges, arising especially from the crossed diagram. It is shown how this allows one 
to recover current single boson exchange potentials in the case where the boson is unique 
(neutral) and couples predominantly in a scalar way with the constituent particles as 
in the Dirac or Klein Gordon equations. A simple explanation for the role played by 
this contribution is proposed. Some attention is given to the massless case, with regard 
to possible m a^ log - corrections to the binding energy, and to the nucleon-nucleon 
interaction case, with regard to the spin and isospin of the exchanged mesons. Some of 



the results presented here may be partly known in various domains of physics, but do not 
seem to be (or have been forgotten) in the domain of intermediate energy physics, at the 
border line of particle and nuclear physics. While they may not be quite new, we believe 
that reminding them here will make the paper more consistent. 



Contribution of the field-theory based single boson 
exchange to the current two-body effective inter- 
action 
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Figure 1: Time-ordered single boson exchange contributions to the two-body interaction with 
indication of the kinematics in the center of mass. 

Usual two-body interaction models suppose the instantaneous propagation of the ex- 
changed bosons, which results in the conservation of the number of particles. In a field 
theory approach, the number of particles is not conserved and the two constituent particles 
may be accompanied by one, two,... bosons in flight (beside the excitation of anti-particles 
ignored in the present work). In this approach, which has a more fundamental basis [Q 
and is much more ambitious, the lowest order interaction of two constituent particles, 
possibly bare, can be obtained by summing the contributions of the two time ordered 
diagrams shown in Fig. 1. In the center of mass system, it reads: 



VEip.Pf) = E ^^ ( r. ,.\. .„ + 



2ujI \E-uol-Ei~Ef E-uol-Ei-Ef) 

In this equation, the summation on the index v is carried over the different species of 
bosons, if any. The quantities g^ and O^ represent the coupling constant relative to the 
interaction of the boson with the constituent particles and its vertex function. This one 
may involve the spin, the isospin and some form factor. It also contains normalization 
factors relative to the constituent particles, . /^ and ./^, so that g will be in what follows 

a dimensionless constant. The quantity, ool, Ei, Ej represent the on-mass shell energy of 
the exchanged boson and the constituent particles. They are defined as : 



Ul = V/i? + fc2. 



E, = Jm? + Pi 



Ef = Jm^ + p}, (2) 



where k, pi, pf represent the S-momenta of the particles under consideration (in the center 
of mass for the two last ones) (see Fig. 1 for the kinematics). As to E, it represents the 
energy of the system we are interested in. Its presence in Eq. (2) is a characteristic 
feature of the coupling of the two-body component to components that are not explicitly 
specified with three or more bodies in the general case. 

When the contribution of Ve to the Born scattering amplitude is considered, the en- 
ergy conservation relation, E = 2E.i = 2Ef, is fulfilled and it can then be checked that 
the two terms in Eq. (1) sum up to reproduce the standard Feynman expression for the 
propagator of a boson. Contrary to this one however, Eq. (1) provides a well determined 
procedure for extrapolating the interaction off-energy shell. The dependence of the in- 
teraction, Ve, on the energy E, which has the same origin as the one present in the full 
Bonn model of the NN interaction |]12| (not to be confused with the other Bonn models 
denoted Q, R, QA, QB, QC, RA, RB |1^, |1^), is at the same time a distinctive feature 
of a field-theory approach, but also a source of difficulties 0. Indeed, inserting Ve in a 
Schrodinger equation leads to solutions that are not orthogonal in the case of states with 
different energies. This violates an important property expected from quantum mechan- 
ics, therefore suggesting that this way of proceeding is wrong or, at least, incomplete. 
While this incompleteness may be remedied by inserting contributions due to the many 
body components that the energy dependence of Ve implies (these may be 2 constituent 
particles plus 1, 2,... bosons in flight), most often, the use of energy dependent potentials 
has been discarded, not without any reason |]l8l . Another way to remedy these difficulties 



is to derive from Ve, Eq. (1), an energy independent, but effective interaction. Many 
works along these lines have been performed in the literature |l^, 20, 21, Q. Here, we 
will proceed in a way which is perhaps more appropriate to emphasize the physics hidden 
behind the above energy dependence HTSf. The interaction Ve, Eq. (1), is expanded in 
powers of l/uj^. Starting from: 

1 1 E-^-^ 

- 2rn 2m +...^ (3) 



2m 2m 



E-ojk-^-^ ^^ "^k 



where the energies of the total system and of the constituent particles are, from now on, 
replaced by the non-relativistic counterparts, one gets: 

Ve = Vo-1{{E-^),V}---. (4) 

2 m 

A different expansion around some energy, Eq, (renormalization point) could also be 
performed. More appropriate to the framework of the renormalization theory, it is less 
convenient here where a comparison with standard descriptions of two-body interactions 
is intended. The difference, which involves higher order terms in Eq. (3), is irrelevant for 
the present work, essentially devoted to a non-relativistic description of two-body systems. 

The first term in the r.h.s. of Eq. (4), which stems from the similar one in Eq. (3), can 
be identified to a usual boson exchange interaction. In momentum space, it reads: 



In configuration space (for tlie simplest case wliere tlie vertex functions 0„ do not depend 
on momenta and wliere form factors are neglected), its expression is given by: 

Vo{r) = -Y.9lOlOl''^, (6) 

The quantity Vi in Eq. (4), which is related to Ve by the relation : 

Vi{k) = — -^, (up to — terms), (7) 



is given, in momentum space, by: 



^iW=EV%^, (8) 

The corresponding expression in configuration space, in the same approximation as for 
Eq. (6), reads: 

Vi{r)-}_^9^0^0^—J ^^^^^ , (9) 

With the expansion of Ve we made together with the assumption of local vertex functions, 
the dimensionless quantity Vi{r) given by Eq. (9) is local. This feature greatly facilitates 
the interpretation of its role but will not be necessarily used throughout this paper. In 
the case of massless bosons, like in QED or in the Wick-Cutkosky model, the integral in 
Eq. (9) diverges logarithmically and some slight modification of the expansion given by 
Eq. (3) should be performed (see Sect. 3.2). 

With the above definitions, the equation that the wave function describing the two-body 
components, ip{r), has to satisfy reads: 

(vo{r) + i{^, (1 + \/i(r))} - E{1 + V^ir))) ^{r) = 0. (10) 

By multiplying on the left by (1 + Vi(r))^^/^ and making the substitution 

^(r) = (l + \/i(r))-i/2^(r), (11) 

one gets an equation that can now be identified to a Schrodinger equation: 

lv{r) + ^-E]<j){r) = 0, (12) 

y m J 

where 

\/(r) = (1 + V,ir)r'/' Voir) (1 + Vir))-'/' + .... (13) 

The dots in Eq. (13) represent a higher order contribution in both Vi(r) and — , whose 
expression, in the limit of a local operator, is given by: 

1 (V\/i(r))2 



4m(l + \/i(r))2' 



Being of the order of relativistic contributions neglected here (Z-diagrams for instance), it 
is not considered in the following. The derived potential, V{r), can thus be approximated 
by: 

V{r) = Voir) - ^{Vi{r), Vo{r)} + ■■■. (14) 

This expansion may be useful for a discussion of corrections order by order. The full 
expression (13) is nevertheless more appropriate in the case where an expansion in terms 
of Vi is not legitimate. 

At the lowest order in Vi, the correction to the potential Vq reads in momentum space: 



^ ^^l + {p^-p)' ^ {fil + {pf-p)2)l)- 

After making the change w ^->- v in the second term and using quantities u]" and uj^r 
defined analogously to Eq. (2): 

< = (/"l + (p. - P)')i ^} = [P^l + {Pf - P)')^ (16) 

the above expression may be written in a simpler form: 

This correction, which contains twice the coupling g"^ appears as a second order boson 
exchange contribution and, for a unique spinless scalar boson, tends to have a repulsive 
character. For consistency, it should be considered with contributions of the same order 
(to be discussed in Sect. 4). The appearance of second order contributions in Eq. (13) may 
look surprising as the starting point given by the field-theory motivated interaction, Eq. 
(1), was a first order one. It originates from the off-energy shell behavior of this interaction, 
which, after elimination, transforms into a higher order correction. This currently occurs 
when changing the degrees of freedom used in describing some system to effective (or 
dressed) ones (see |^3|] for another example in relation with the NN interaction). Notice 
that Vo(r) and Vi{r) in Eq. (14), or factors O^ O^ and Ol O^ in Eq. (15), may not always 
commute. This respectively occurs for momentum dependent (non local) interactions or 
spin-isospin dependent ones. In the following, we will sometimes disregard these cases, 
essentially for simplifying the discussion and providing a more sensible presentation of the 
role of corrections associated with Vi{r). 

The operator, Vi{r) {Vi{r)/{1 + Vi(r)) more exactly) has a simple interpretation. It 
represents the probability that the system be in a component where the two constituent 
particles are accompanied by a boson in flight. The interaction V{r) in Eq. (12) therefore 
corresponds to the interaction due to the component with two constituent particles only. 
Consistently with contributions of diagrams of Fig. 1 retained until now, it does not 
contain any interaction of these particles while a boson is in flight (see Fig. 2 for a 
graphical representation). As to the reference to dressed particles, it makes it clear that 
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Figure 2: Time-ordered boson excliange contributions to the two body interaction, including 
some interaction between the two constituents (symbolized by the oval). 



higher order corrections should appear in the effective potential, as mentioned above. 
These corrections have some relationship to the renormalization presented in a more 

Let's also mention the consistency of: 
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S|. Let's also mention the consistency of: i) the 
above interpretation of Vi{r), ii) the general definition of the normalization for an energy 
dependent potential like that given by Eq. (10) and iii) the relation given by Eq. (11): 



df'^\r){l - —^) = / dfilj'{r){l + Vi{r)) = / df(f)\r), (18) 

Not surprisingly, there is a relationship to the renormalization in the field of the many 
problem. The quantity, (1 + Vi{r))~^, is quite similar to the spatially dependent spectro- 
scopic factor, which is introduced in the domain of mean field approximations to account 
for the effect of degrees of freedom (nuclear correlations) that have been eliminated (see 
ref. 1^^ for a review). Its effect, or what accounts for it to be more correct, turned out 
to be essential to preserve the orthogonality of many body wave functions when incorpo- 
rating contributions due to the same degrees of freedom [^ . 



At this point, it thus appears that the interaction V{r), which is obtained from a 
field-theory based single boson exchange interaction does not identify to the single boson 
exchange potential Vo(r) which, most often, is referred to and is used in the Schrodinger 
equation, or an equivalent one. 

In the case of the strong A^A^ interaction, the above problem may not be a serious source 
of concern. Its short range part has some phenomenological character and through the fit 
of the potential parameters to experimental data, it may be possible to account in some 
part for the correction to Vo(r) approximately given by |{Vi(r), Vo(r)}, which also has a 
short range character. Only the longest range part of this term which corresponds to two 
pion exchange and has a complicated spin isospin structure, may not be well reproduced. 

In the case of the Coulomb potential, where the fundamental interaction is known, there 
is no uncertainty similar to that affecting the short range part of the A^A^ interaction. As 
contributions provided by the term, |{Vi(r), Vo(r)}, are known to be absent from the 
energy spectrum of atomic bound states involving electrons or other charged particles, 
the derivation of the potential given in this section has to be completed in any case. This 
will be done in Sect. 4, while contributions will be estimated in the next section. 

It is worthwhile to mention that removing the energy dependence of the potential Ve is 



quite similar to the Foldy-Wouthuysen transformation for the Dirac equation. Like this 
one, it is not merely a mathematical operation. From a state which is a superposition of 
a two-body and a two-body plus one boson in flight components, which respectively play 
the role of the large and small components in the Dirac spinors, it allows one to build a 
state with a unique two-body (dressed) component. This one represents a frozen coherent 
superposition of the bare components, with the consequence that the constituent particles 
have now acquired in the transformation an effective character. The change in the nature 
of the constituents is better seen by working in Fock space |]2S[ (see also the appendix for 
an example). Quite similarly to the Foldy-Wouthuysen transformation, the developments 
presented in the present section make sense around some given energy (states close to 
threshold in the present work). On the one hand, the choice of this energy, which plays 
the role of the renormalization point in the renormalization theory, could be optimized. 
On the other hand, by removing higher order energy dependent terms in the expansion 
given by Eq. (3), one can enlarge the energy domain at will, as far as one does not cross 
the inelasticity threshold where the frozen degrees of freedom may show up. While doing 
so however, the effective (or renormalized) interaction V{r) in Eq. (13) acquires a non 
local character whose role will steadily increases when approaching this threshold. 

The contribution considered here may have been considered as a relativistic effect and, 
as such, as a new contribution to be accounted for. It implies a non-instantaneous prop- 
agation of the exchanged boson, but this feature is also pertinent to a two-body system 
that, under the effect of some interaction, undergoes a transition from a low momentum 
component to a high momentum one and back, or a two constituents system (NN for 
instance) that undergoes a transition to a component involving their internal excitation 
(AA) and back. These processes, which can give rise to an effective interaction, are typ- 
ically part of a non-relativistic description of a two-body system. Contrary to current 
relativistic effects that are of the order (f /c)^, they have a static character and don't 
vanish when the mass of the constituents goes to infinity. 

The current relativistic effects and the contribution of interest here can be related to 
different contributions in the expression which, for instance, allows one to calculate some 
amplitude from an interaction acting on an input amplitude: 

'^'^J "^^^ ((|Po + qoV -m^- ....)((|Po - qoY - m^ - ....) ^7^' ^^^^ 

In this expression, where one recognizes the propagators of the constituent particles as well 
as the propagator of the exchanged boson, go and Pq represent the time component of the 
4-momentum carried by the boson and the system under consideration. Contributions to 
the integral come from the poles in the integrand. The poles arising from the constituent 
propagators generally imply a small value of go? determined by the energy transfer of 
a constituent in the initial state to a constituent in the final state. They give rise to 
the above mentioned relativistic effects of the order o;"'^. The pole due to the boson 
propagator implies a value of go that by no means is small. It is responsible for the 
effect mainly considered here, of the order a;~^. While it is conceivable that the first effect 
disappears in the limit of an instantaneous propagation of the boson, which supposes that 
go can be neglected in the boson propagator in Eq. (19), this is difficult to imagine for the 
second one. It is also noticed that the two contributions are on a quite different footing. 
The second one disappears when the boson propagator in Eq. (19) is replaced by the 

9 



propagator (/x^ — g^ + p-i )■, which is equivalent for a physical process (Born amplitude). 
On the other hand, as it can be seen from the expression of the denominator in Eq. (3), 
the effect which is considered here is of the order /x/m, i.e. a recoil effect. 



3 Estimates of the contributions due to the term 

\{VuV,} 

We present in this section results of perturbative calculations of two-body wave functions 
with and without including in the Schrodinger equation the contribution due to the term 
|{^i;K)}- At the same time, we will discuss the relationship of the wave functions so 
obtained with those calculated from an energy dependent interaction. The aim is to show 
that a large part of the physics which is at the origin of this energy dependence can be 
accounted for in an energy independent scheme (where the underlying degrees of freedom 
have acquired an effective character). In practice, this will be done for the NN interaction 
case and the Wick-Cutkosky model. 

The wave functions of interest are given in momentum space by: 

M?^ = -^ / (^ ''»(?. P') ^o(p'). (20) 

where Vip^p ) differs from Vo(p, p ) by the correction given by Eq. (15). The quantity /t^ 

2 

is related to the binding energy by the relation \E\ = —. By comparing the results, we can 
estimate the effect of this correction. In principle, if 0o(p) is chosen as the solution of the 
Schrodinger equation with the potential Vq{p, p'), (poiP) is equal to 0o(p) (strictly speaking, 
this is not anymore a perturbative calculation but, rather, a consistency check). In some 
cases however, some terms in the original potential model that was used to calculate (po{p), 
may be neglected. It is then more appropriate to look at the effect of the correction by 
comparing 0(p) and 0o(p), rather than (f){p) and (po{p), hence the introduction of (po{p)- 
On the other hand, while a schematic description of effects under consideration here has 
been given in the previous section, in practice some of the calculations presented in the 
present section, with the Bonn models, have been performed with the full complexity of 
the interaction. This involves in particular the non-locality due to describing spin 1/2 
particles by Dirac spinors and normalization factors appropriate to the definition of the 
potential V (together with the corresponding equation). For these cases, the numerical 
accuracy has been pushed to the point where one can reproduce the numerical deuteron 
wave function at a level where the discrepancy cannot be seen in the figures presented 
below. 

Other quantities of interest correspond to perturbative solutions of the Schrodinger 
equation where the potential Vo(p,p') is replaced by the energy dependent potential given 
by Eq. (1): 

m = —7^ / 1$^. Ve{p,p') Up)- (22) 



p'^ + K^J (27r)3 
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Making an expansion up to the first order term in Vi, the above equation may also read: 
~ m f dp 

{vo{p,p) + i(^!±^ v,{p,p') + vMv) ^^^±^)\ Mp') 

\ 2 m m J 

^ ^ '^^' V,{pJ')Up)- (23) 



2 7 (27r)3 

Equation (22) (together with Eq. (23)) differs from Eq. (21) by the fact that the first 
one involves the full energy dependent interaction V^, while the second one involves the 
corresponding energy independent potential. It may be of interest to compare the outputs 
of these equations. It is remembered that 0o(p) and ^o(p) differ by an operatorial factor, 
(1 + Ki)2, see Eq. (11). Two cases are to be considered: 

i) The first one assumes that ipo (and not 0o) is solution of the Schrodinger equation with 
the potential Vq. One thus gets: 



This is formally close to the choice made in refs. |TI], |^. 

ii) The second case assumes that 0o(p) is solution of the Schrodinger equation with the 
potential Vq. According to Eq. (11), the zeroth order wave function to be inserted in Eq. 
(23) is related to the above one by the relation %pQ = {1 + V^i)~20q. The wave function, 
Tp{p), then reads: 

;yi{p,p) ipoip') 



2 J (27r)= 
^(l-ivi)0(p)- (25) 

This relation is in agreement with Eq. (11), which relates the wave functions of the energy 
dependent and energy independent schemes, ip ~ (1 + Vi)~20 ~ (1 — ^Vi) (p. 

We now consider successively the Nucleon Nucleon interaction and the Wick Cutkosky 
(or QED) models. 



3.1 Nucleon Nucleon case 

We estimated corrections to the wave function arising from the term |{Vo, Vi}, starting 
from a zeroth order wave function given by the Paris model [0. This was done per- 
turbatively, using expressions given in the introduction of this section. In a first step, 
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Vq was approximated by a sum of a few well known meson exchanges {ir, p,uj,a), whose 
contributions in the standard non-relativistic limit were retained. The various parame- 
ters were chosen so that to reproduce the deuteron wave function of the model at large 
distances and as accurately as possible at short distances. Results obtained in this way 
have shown that the deuteron S wave function was strongly reduced, more than obtained 
in refs. [|TU|, 0|. The D wave was less changed, resulting in a relative enhancement of the 
D wave with respect to the S wave, qualitatively in agreement with the findings of these 
works. Due to the large effects, and to avoid to be misled by a hidden bias, we turned to 
a more complete calculation, with Vq taken as identical to the potential used to calculate 
the deuteron wave function. In this order, we use meson exchange models which, contrary 
to the Paris model, allow one to easily derive the operator, Vi. The choice of the Bonn 
Q model will make possible a comparison with earlier work JlOl. As this model provides 



a wrong prediction of the mixing angle, ei, for energies beyond 100 MeV, we also use the 
Bonn QB model which provides a correct prediction with this respect and is generally 
considered as a quite reasonable model. 

In Fig. 3, we first show the effect of introducing somewhat arbitrarily extra normaliza- 
tion factors J^ J^ in the integrand of Eq. (20). Calculations are performed with the 
Bonn QB model for both the deuteron S and D waves. This can provide an estimate of 
the size of typical relativistic effects. The largest part, 20% in the range p = 1 GeV/c 
rT| , is given by the factor, J^-, which can be factorized out in Eq. (20). As the pres- 



ence or absence of such a factor is dictated by unitarity conditions and is accounted for 
in many nucleon nucleon interaction models, the effect may not be a real one however. 



Some effect is also due to the other factor, w^, whose relevance was mentioned in |I4 . 
It gives rise to a decrease of the S wave of 4% at small momenta, which has the size 
of corrections generally expected from relativistic effects involving the internal structure, 
and of the order of 10-20% in the range p = 1 GeV/c. The effect of the factor, (1 -|- Vi)~2 , 
relating the wave functions in the energy dependent and energy independent schemes, is 
also shown in Fig. 3. Due to the difficulty in performing a full calculation, which has to 
do with the non-locality of the operator, we can only display the effect due to first and 
second order corrections in Vi. This however has the advantage of providing an estimate 
of some contributions neglected here. As expected, wave functions at small momenta are 
not modified. At intermediate momenta, some decrease is observed for both the S and D 
waves, which is compensated at the level of the norm by the contribution due to the me- 
son in fiight content of the order of 3-4%. The results roughly agree with those obtained 
in ref. [^, which were calculated exactly but with a localized form of Vi, or with that 
obtained for the full Bonn potential [0. Altogether, results presented in Fig. 3 show 
that neither relativistic effects, nor those dealing with the relation of the wave functions 
of the energy dependent and energy independent schemes can explain a large change of 
the deuteron S and D waves at small momenta (and equivalently at large distances). 

In Figs. 4 and 5, we show the results of a perturbative calculation tending to estimate 
the effect of corrections to the potential due to the extra term, \{Vq, Vi\ (dotted line). In 
this order, we insert as inputs in Eq. (21) the wave functions and potentials corresponding 
to the Bonn-Q and Bonn-QB NN interaction models. The resulting wave functions for 
these two models and for both S and D waves are drawn in Figs. 4 and 5, together with the 
unperturbed ones. It is observed that the global structure of the wave functions remains 
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Figure 3: Sensitivity to relativistic kinematics and underlying schemes (dependence or inde- 
pendence of the total energy). The input wave function (continuous line) and the potential are 
those of the Bonn-QB model. The S and D components are presented respectively in the top 
and bottom parts of the figure. To better emphasize the small momentum behavior, we provide 
in the insets the ratio of the various calculations to the exact ones. 

essentially unchanged for the Bonn-Q model and the S wave of the Bonn-QB model. 
The position of the minima is only slightly shifted. Quantitatively, the most important 
feature concerns the decrease of the S wave at low momenta. This one is largely due 
to a decrease of the contribution produced by the tensor force acting on the deuteron 
D wave. The effect, which is better seen in the insets of Figs. 4 and 5, is given by a 
factor 0.73 for the Bonn-Q model and 0.63 for the Bonn-QB model. It shows that the 
wave function accounting for the renormalization of the interaction by the implicit meson 
in flight probability, see Eq. (13), strongly differs from the original one in a momentum 
domain, p ~ k (= 0.232 fm~ ), which is essential for the deuteron description. In all 
presented calculations, the D wave at small momenta is much less affected, resulting in 
an enhancement of this one when the complete perturbed wave function, Eq. (21), is 
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Figure 4: Effect of the ^{Vi,Vo} term, Eq. (14), on the calculation of the S component of 
the deuteron wave function. The input wave function (continuous line) and the potential are 
those of the Bonn-Q and Bonn-QB models (respectively upper and lower parts of the figure). 
The dotted and dash-dotted lines correspond to outputs of the energy independent and energy 
dependent schemes respectively. The low momentum results are shown in the insets as a ratio 
to the exact results. 



normalized to unity. The comparison of the results obtained with the two Bonn models 
does not show significant qualitative differences at low momenta, despite the fact that the 
model Bonn-Q does badly for observables in relation with the tensor force. This rough 
agreement indicates the degree of reliability one can give to the estimate of the effects 
under consideration in the present work. 



The comparison of present results with those of refs. ||T0|, evidences some difference: 
their suppression of the S wave, roughly a factor 0.8, is not as large as here. Furthermore, 
the "relativistic" wave functions in the first reference, /i and /2, significantly differ from 
their non relativistic limits respectively given by the usual S and D waves: a new zero 
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Figure 5: Same as in Fig. 4, but for the D wave. 



appears in /i while none is seen in /2 in the range extending top = 1.5 GeV/c (~ 7.5 fm" ). 
To identify the origin of these differences, which may have many sources and among them 
genuine relativistic effects not accounted for here, we did a calculation whose spirit is 
closer to theirs. Using Eq. (22), we took the energy dependent potential given by Eq. (1) 
together with the unperturbed wave function given by the Bonn-Q model. Calculations 
corresponding to an expansion of Ve up to the first order term in Vi have also been 
performed. They may allow one to test the validity of this approximation, which was 
made in deriving the potential appropriate to an energy independent picture, Eq. (13). 
Results are presented in Fig. 6 for both the S and D waves, together with the unperturbed 
wave functions. As for calculations presented in Figs. 4 and 5 for the same Bonn model, 
we do not observe major changes in the global structure of the wave function. The position 
of the zeros is only shifted a little more, as is the position of the same zeros in Fig. 3 
for the Bonn-QB model. Quantitatively, the S wave at low momenta is also reduced, by 
a factor 0.89 for the full operator Ve and a factor 0.85 for its expanded form limited to 
the first order term in Vi. The decrease compares to that obtained in refs. |jlO|, 0, but is 
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Figure 6: Effect of the energy dependent part of the potential, Eq. (1), on the calculation of 
the S and D components of the deuteron wave function with various approximations: full kernel 
(dotted line, Ve) and limited expansion up to the first order in Vi (dash-dotted line, L.E.). The 
input wave function (continuous line) and the potential are those of the Bonn-Q model. The S 
and D components are presented respectively in the upper and lower parts of the figure. The 
low momentum part is shown in the insets as a ratio to the exact results. 

slightly smaller. We checked that this can be ascribed for a large part to the kernel of the 
interaction in the underlying light-front approach, which has some similarity with ours, 
Eq. (1), but also some difference. In any case, the decrease is significantly less than the 
one obtained previously in the energy dependent scheme (0.73), suggesting a bias in one 
of the approaches whose results are compared. 

The difference in the reduction can be traced back to the extra contribution due to the 
term | VqVi, which is present in Eq. (21), but absent in Eq. (24). This term is a priori 
different from the other one, | Vi Vq, which is present in both equations. They become 
equal to each other in the approximation where Vq and Vi are local and commute. This 
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can explain the rough factor 2 difference between corrections calculated from Eqs. (21) 
and (24). As explained in the introduction to this chapter, it would be more appropriate 
in a calculation of the wave function for the energy dependent scheme to use the function, 
■?/;o — (1 + Vi)~20Qj as an unperturbed one. When this is done, the above discrepancy 
disappears, the remaining difference being provided by the operator, (1 + Vi)~2 ~ 1 — |Vi. 
Results calculated in this way at the first order in Vi are presented in Figs. 4 and 5 (dash- 
dotted line). As expected, the wave functions calculated perturbatively in both schemes 
coincide at small momenta. The differences at higher momenta are similar to those 
already mentioned in relation with the discussion of results presented in Fig. 3 (effects of 
the factor (1 + Vi)-^ ~ 1 - i^i). 

None of our results for the Bonn-Q model presents a qualitative change of the wave 
function similar to that obtained in ref . [|ry] , namely an appearance and a disappearance 
of a zero in the S and D waves respectively (below p = 1.5 GeV/c ^ 7.5 fm~ ). As 
mentioned by the authors, this feature has been considered as a genuine feature of the 
relativistic wave functions in their approach. Surprised by the fact that we could reproduce 
other features of their calculation, especially the decrease of the S wave at low momenta, 
translating into a relative increase of the D wave after normalization, we performed further 
checks. We found that the parametrized wave function given in ref. [jl2|, whose behavior 



at high momenta differs from what is theoretically expected, was not accurate enough 
to allow the authors to make a so strong statement. The wave functions that come out 
from a consistency check, Eq. (20), differ from the input by an amount which, in the 
range p > 0.8 GeV/c (~ 4.0 fm~ ), is of the same order as the effect they were interested 
in. Furthermore, the potential effectively used differs from the original Bonn-Q model in 
the treatment of the p meson nucleon tensor coupling. This involves a term containing 
the time component of the 4- momentum carried by the meson, go- We re-calculated what 
should be the reference wave function for their study, using the parametrized wave function 
on the one hand and their p meson nucleon coupling. Results, which are presented in Fig. 
7 (dash-dotted line), are qualitatively quite close to their so called relativistic components, 
/i and /2. The extra zero in the former and the absence of a zero in the latter, in the 
range p < 1.5 GeV/c (~ 7.5 fm~^), thus appear as the consequence of approximations in 
the calculations. Other details relative to the role of various ingredients are shown in the 
figure. 



3.2 Wick-Cutkosky model, QED 

The exchange of zero mass bosons, which applies to the Wick-Cutkosky model or to QED, 
especially the Coulomb potential, deserves some attention. Indeed, the quantity, Vi{r), 
defined by Eq. (9) is then scale independent and, being proportional to J dx{smx) /x^, 
is logarithmically divergent, requiring to come back on the definition of Vi{r). 

Two questions are of interest here: the presence of a log - corrections to the binding 
energy and some correction of the order p/m to the wave function [Q which, both, 
have been attributed to a relativistic approach. 

To determine Vi(r), we started from an l/uk expansion of the propagator {E — ojk — 
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Figure 7: Consistency clieck and sensitivity to the deuteron input wave function for the Bonn-Q 
model (upper and lower parts for the S and D waves respectively). The continuous and dotted 
lines respectively represent the exact (numerical) and parametrized wave functions (denoted 
exact and param. in the figure). The wave function recalculated using the exact wave func- 
tion cannot be distinguished from the original one while that one using the parametrized wave 
function (param.-recalc), represented by the short-dashed line, evidences some discrepancy. The 
same wave functions recalculated using the p meson nucleon tensor coupling employed in ref. [|l^] 
are respectively represented by the long-dash and dash-dotted lines (denoted exact recalc.(with 
qO) and param.-recalc. (with qO)). 



2m ~ 2m)~^' ^^- ('^)' "^ith the idea to only retain the first order term in {1/ujkf'. The 
presence of a singularity in the higher order terms when the boson mass is zero, which 
implies that tUfc — ;► for A; — > 0, supposes to proceed differently. The difficulty can be 
overcome by removing from the full propagator, Eq. (3), that part which gives rise to the 
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usual zero mass single boson exchange contribution to the two-body interaction: 

1 1 E 



2m 2m 



E-oJu-i-'i '^^ uu{uu-E + i-+'i 



(26) 



Now, the extra quantity at the denominator of the second term in the r.h.s. of the 
above equation, —E + §^ + 2m^ '^^^ ^^ bounded from below by the absolute value of the 
binding energy of the system under consideration, ■^, where k = Jm\Eh^(,]. This, which 
corresponds to summing up an infinite set of terms in Eq. (3), makes the second term in 
the above equation less singular than in the former, allowing one to deal with it rather 
safely. Being positive, the extra terms, ^, don't introduce poles and, on the other hand, 

2 

their size is expected to be the same as that of the term — , which we will account for by 
multiplying this one by a factor 2. After making these approximations, the developments 
from Eq. (3) to Eq. (9) can be performed again with the result: 

2 f dkeUkr) 



7T J ^i(^j^ + 2^y 

2^ r dkjojkr) 



which is seen to be a function of the dimensionless variable 2— r. The coupling has 
been particularized to the case of a spinless particle (the time component for the photon) 
and the coupling g^ has been expressed in terms of the coupling, a, often used in QED, 
g"^ = Alia. In terms of this coupling, the quantity k for the ground state of the system of 
two equal mass particles is given by /t = ^. 

The small distance behavior of Vi(r) is dominated by the logarithmic divergence of the 
integral in Eq. (27): 

Vi(r).^o = - a (l - 7 + log(^)) + 0(— ). (28) 

where 7 represents the Euler constant, 7 = 0.577... At very large distances (r >> ^), 

Vi{r) behaves like: 

in 
Viir)r^oo = a-^. (29) 

Zk,^ r 

A useful and approximate interpolating expression may be given by: 

V^ir) c^lalog(l + TT^e^"^) • (30) 

TV \ Zk"^ r / 

For small r, it is exact and for r -^ 00, it provides a slight underestimate by a factor 
ei-^/(7r/2) = 0.97... 

Solving the Schrodinger equation with the potential V{r) together with Vi{r) given by 

Eq. (28) may be done approximately by noticing that Vi{r) is a smoothly varying function 

in the range r ~ ^ >> — , corresponding to a << 1. It can then be approximated by a 

constant number: 

2 1 

Vi{r)rc^^-i = -alog-. (31) 

n a 
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a 



With the above approximation, the correction to the potential given by |{Vi(r), Vo(r)} can 
be incorporated in the Coulomb potential, Vo(r) = — -, whose coupling is renormalized 
as follows: 

a(l alog-V (32) 

The correction to the binding energy at the lowest order immediately follows: 

a^ a^ / 4 1\ 

Eb.e. = -m —^ -^ Eb.e. = -m — — 1 a log - . (33) 

In all the previous discussion, we have been mainly concerned with the largest correc- 
tion, of the order a log -. We have not paid much attention to the precise value of some 

factors. Quantities such as — or r have been replaced by approximate values in getting 
Eqs. (31, 32), which may affect the coefficient of ^ in the log term. This can be seen to 
correspond to a correction to the binding energy of the relative order a in Eq. (33), but 
does not change the correction of the order a log -, which dominates in the limit a — >■ 0. 
Going beyond requires some care, but one can anticipate that the replacement of the 
factor 2— appearing in Eqs. (27-30) by the potential a/r, would provide a useful starting 
point for a discussion. With this approximation, motivated by the origin of the factor, 
—E + —, which, acting on the wave function, provides the potential a/r, V\{r) becomes 
a constant. Expressions with increasing accuracy would be successively given by: 

2 1 

Vx{r) ~ - a log - [up to 20% for a < 0.1], 
n a 

2 1 

~ - a (1 - 7 + log -) [up to 25% for a < 0.5], 
vr a 

2 / e^-'^\ 

- - a log 1 + [up to 4% V a]. (34) 

vr Y a J 

The first expression is identical to Eq. (31), which was obtained on a slightly different 
basis. The second and third ones could be used to improve Eqs. (32) and (33), when 
the first one is out of its range of validity (limited to a < 0.1). A further improvement 
consists in expressing the quantity, 2—, in terms of the renormalized potential, aeff/r. 
The effective coupling, ag//, which should replace the factor a in the log term appearing 
in Eq. (34), is determined by the equation: 

"e// = Y^^ = «(l-"e//'^)' (35) 

where J = - log -^—, - (1 — 7 + log -^—), - log(l + - — -) or a better expression depending 

^ J J ^J J ^J J 

on the desired accuracy. 



The correction to the binding energy in Eq. (33) is identical to that found in [0. 
This result, which has been obtained here in a picture that is not especially a relativistic 
one, tends to show that the correction has essentially its source in the field theory which 
underlies the present approach as well as that followed by Feldman et al. ^ and many 



other authors [15, M]. Results for the binding energy of the lowest states calculated with 



the renormalized interaction have also been obtained |3T| . They explain the bulk of the 
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departures from the instantaneous approximation results derived either from the Bethe- 



Salpeter equation [g2| or from energy dependent approaches |^, ^ 



It is also interesting to consider various wave functions corresponding to different ap- 
proximations or different schemes. In comparison with the NN interaction case discussed 
previously, one can go a step further here by using the renormalized interaction which, 
in a first approximation, differs from the bare one by changing the coupling a into Oe//, 
Eq. (35). In Fig. 8, we present wave functions pertinent to the lowest s-state. These are 
the usual non-relativistic wave function: 

5 

(f)o{p) = V^ , 2 , °9\9 ^i^h '^0 = 77"^ a, (36) 

the non-relativistic wave function with the renormalized interaction: 

(j){p) = v^^-^— -^ with K = -m aeff, (37) 

the wave function corresponding to the energy dependent scheme obtained by using Eq. 

(ll),^o = (l + ^i)-^0o: 

Mp) = J^^m, (38) 

V a 

and the same wave function calculated perturbatively using Eq. (22), of which analytic 

expression is given by: 

. 5 , 

Ak2 fa, 



"^^^ = ^WTv^- i 



eff 
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a 



K, + p^)"^ V a 
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1 I P. I «^^+P^ ^'^ /o K2+p2 

-^ ^ m ^ 4m2 V ^ ^2 
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(39) 



where the square root (and log) functions should be appropriately changed when their 
argument gets negative (or imaginary). 

In writing Eq. (39), we have emphasized at the first line factors which together form 
the unperturbed wave function, ipoi^p), see Eqs. (37) and (38). The factor that accounts 
for the improvement obtained by using the complete kernel of the interaction, Eq. (1), is 
given at the following lines. For the numerical study, whose results are presented in Fig. 
8, the value of a has been taken to be 0.5, which corresponds to ae//= 0.32. 

Not surprisingly, 0o(p) and 0(p) differ significantly, showing the indirect importance of 
the energy dependence of the interaction, Ve, which was accounted for by employing a 
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Figure 8: Wave functions for the Wick-Cutkosky model in different approximations: from the 
bare Coulomb potential, Eq. (36) (dash-dotted line); from the renormaiized potential, Eq. (37) 
(continuous line); zeroth order for the energy dependent scheme, Eq. (38) (dotted line); and the 
corresponding one from a first order perturbative calculation, Eq. (39) (dashed line, essentially 
on top of the previous one)). Units are expressed in terms of the constituent mass, m. The 
numerical values of the wave functions, normalized as in Eqs. (36, 37), have been divided by 10. 

renormaiized interaction in the energy independent scheme. The wave function appropri- 
ate to the energy dependent scheme, '?/'(p), is quite close to the corresponding unperturbed 
one, ipo{p}- The discrepancy, which can hardly be seen in the figure, is less than 1% in 
the domain of momenta covering 95% of its contribution to the normalization. This in- 
dicates that the way ttg// has been derived is essentially correct. These wave functions 
strongly differ from that one obtained with the bare potential in the energy independent 
scheme, (poip)- A similar result was obtained in the light front approach in refs. [|, |^ 
and most probably with the trial wave function employed in ref. [^|. The point is of 



relevance for works that have some relationship to the present one [^], where only an 
effect in the high momentum range was emphasized. This second effect, which implies 
corrections of the order p/m, is also present here (see below). It certainly explains some 
discrepancy between ip{p) and its unperturbed counterpart ipoip), but the effect shows up 
at high momenta, outside the range shown in Fig. 8 (4% at — = 0.5, 10% at — = 0.8 
and 20% at — = 1.3). We also looked at the contribution that the operator Vi provides 
for the normalization, see Eq. (18). The operator has been chosen according to calcula- 
tions presented in this section. It is given by Eq. (27) together with the full propagator 
appearing in Eq. (26) (instead of —-q§-i which it is equal to in some limit). Discarding 
any theoretical improvement, the contribution was found to be close to what was ap- 
proximately expected, 1 — -^^ (0.345 instead of 0.360). This indicates, in a particular 
case, that the derivation of the effective interaction, V{r), Eq. (13), its solution, 0(r), 
Eq. (12), the relation of this solution with that for the energy dependent scheme, ip{r)., 
Eq. (11), and the normalization condition, Eq. (18) are essentially consistent with each 
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other. The remaining discrepancies, which are at the level of a few % in comparison with 
the major effect of 30-40% related to the renormalization of the interaction, are ascribed 
to the approximations made in dealing with the denominator of the second term in the 
r.h.s. ofEq. (26). 

We now consider more closely how the wave function of the energy dependent scheme 
given by Eqs. (1, 10) relates to the energy independent one given by Eqs. (12, 14) in the 
high momentum range. This relation in the configuration space, 

^(f) = (1 + \/i(r))-i/2 0(f) ~ (1 - i \/i(r)) 0(f), (40) 

reads in momentum space: 

ij{p) ^ <p{p) - 1 I ^ V,{p,p') <l>{p') ^ <j>{p} + 5m, (41) 

The high momentum regime of this equation is obviously that one where relativistic effects 
are expected to show up. 

Neglecting a possible renormalization of the interaction as given by Eq. (13), or as- 
suming it can be accounted for effectively using Eq. (34), the solution of the Schrodinger 
equation for a Coulomb or Yukawa type interaction, 0(p), obeys the equation: 

,, . ATram r dp' 1 ,,^,, ,,„, 

p^ + K,^ J [Ztt)'^ [P — P ) + /i 



while Vi{p,p') is given by: 



^i(P>"") = T— 4^^-^. (43) 

{{p- P ) + fJ' )^ 



In the limit of large p, Eq. (42) indicates that (f){p) behaves like: 



, / -A 47rQ; m f dp' , , ^,, , ^ ^ , 

^(P).-~ = ^3-ip^^(p'), (44) 



where, up to some factor, the integral represents the configuration space wave function at 
the origin. The second term in Eq. (41) can be similarly calculated, using the expression 
(43)ofV^i(p,p'): 

1 Ana f dp' 

2^ J (2^ 

Gathering results of Eqs. (44) and (45), one obtains: 



^^(p),-c„=-i^ 7^^(P'). (45) 



i^{p)p^oo ^ (1 - -— ) 0(p)p^oo, (46) 

2m 

This result supposes that Vi is small enough so that the first order expansion given in 
Eq. (40) is valid. It is actually limited to the range k < p < m and neglects higher order 
effects in the coupling, a. A result which is better with both respects, but limited to a 
zero mass boson, is given by Eq. (39) (see the second line). 
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The above effect may be compared to that one attributed to relativity in |TH]. In 
this work, the hght front wave function was found to be expressed in terms of the non- 
relativistic one as: 

Ep 

where n represents the orientation of the hght front. In the non-relativistic limit, where 
— < 1, and after making an average over the orientation of n, the above result reads: 

Hp)i.f. ^ (1 - £^) '^(P")--' (48) 

which is identical to what was obtained above by an approach of which full relativity 
is absent. This indicates that the correction of the order ^ is not characteristic of a 

m 

relativistic approach, but rather of the field theory which underlies Eqs. (1, 10, 11) on the 
one hand and (46) on the other. Relativity in this last equation would imply higher order 
terms in — , for instance in the replacement of ??i by Ep or in the non-isotropic dependence 
of ip{p)i.f. with respect to the orientation of n. In both cases, they will manifest by effects 

2 

of the order ^, as usually expected. Notice that the above — correction has for some 
part an elusive character. In the low momentum domain, it disappears from a correct 
definition of the norm, see Eq. (18), in the same way that renormalization effects do 
not show up in the calculation of various observables at the lowest order in 1/u;^ p9[. It 
may nevertheless show up in the ultra relativistic regime where — >> 1 |^ but, most 
probably, this is a feature pertinent to the energy dependent scheme. 



4 Extra contributions due to two boson exchanges 

We consider in this section the role of genuine two boson exchanges with respect to the 
correction brought about by the term in Eq. (14). We derive first the general expression 
for these contributions and subsequently look at particular cases depending on whether 
bosons are neutral, charged or massless. 

4.1 General case 

The contributions, which have the same order as the correction given by Eq. (15), and 
are due to two boson exchanges, are calculated using standard perturbation theory. Only 
those which have some relevance for our purpose are retained. The corresponding time- 
ordered diagrams are shown in Figs. 9, 10, 11. They involve a common normalization 
factor, l/{uiUJf) , a common factor corresponding to two boson propagation, l/{uJi + ujf)), 
and two factors corresponding to single boson propagators, l/{uji) or/and l/{ujf). 

The first ones. Fig. 9, have a box type character, and have a topological structure 
similar to that implied by Eqs. (15, 17): 
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Figure 9: Selected contributions to the two-body interaction due to two boson exchange (box- 
type diagrams). Apart from the common factor 1 / (ujiu; f{uJi + to f)) , they involve a factor l/uJitOf. 

The second ones, Fig. 10, have a crossed type character. They involve the two time- 
ordered diagrams corresponding to Fig. 1, while another boson is exchanged. They 
provide most of the contribution to the interaction of the two constituent particles while 
a boson is in-fiight. Fig. 2. Their contribution is given by: 

The third ones. Fig. 11, have also a crossed type character. Except for a different 
spin-isospin structure, their contribution is identical to that shown in Fig. 9 and is given 
by: 

Notice that the ordering of Ol, O^, O^ and O^ in Eqs. (50, 51) is different from that in 
Eqs. (15, 49). 

In deriving Eqs. (49-51), we neglected contributions involving the energies of the con- 
stituents, consistently with our intent to only consider second order effects in boson ex- 
changes. They may obviously be accounted for, which supposes to extend the analysis we 
developped for one boson exchange, Eq. (1), to two boson exchanges. 

On the other hand, we did not consider the excitations of a constituent and its anti- 
particle (negative energy state). These ones, which imply Z-type diagrams (see Fig. 12), 
should be in any case accounted for, but having a different mathematical structure (the 
energy of the excited pair at the denominator and different vertex functions), they don't 
interfer with the discussion presented here. In the case of infinitely massive constituent 
particles and exchanged bosons coupling in a scalar way to them, their contribution should 
vanish. 

The momenta in Eqs. (49-51) refer to those of particle 1, thus introducing an apparent 
asymmetry between particles 1 and 2. It can be verified that for the first contribution, 
AV7, the reference to particle 2 would be accounted for by the change of sign of the 
momenta, Pi, pf and p, which does not provide any change of these equations. For the 
contributions produced by the crossed diagrams, AV// and AVjn, Figs. 10, 11, the 
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Figure 10: Selected contributions to the two-body interaction due to two boson exchange 
(crossed-type diagrams). Apart from a common factor l/{uJiUJf{uJi +0;/)), they contain a factor 
l/ojf (a, b) or l/cj? (c,d). These provide the dominant contribution to diagrams displayed in 
Fig. 2. 

argument is not so straightforward, due to the fact that the momentum of the second 
particle in the intermediate state is not -p, but p — Pi — Pf (in the cm. system). The 
symmetry is recovered by making a change of variable. 



4.2 Neutral and spinless boson case 

The results for the case of the exchange of a neutral boson coupling in a scalar way 
and non-relativistically to constituents stem from Eqs. (49-51). This applies especially 
to the a and u mesons coupling to nucleons and to the photon coupling to charged 
particles (Coulomb part). In both cases, their contribution provides a large, if not the 
dominant part of the interaction. Taking into account that the vertex functions, O^ and 
Ow, commute, it is easy to sum up the different contributions, AV[, AVu and AVm, with 



the result that the common factor, 



[ujr 



+ ^}) 



AVj + AVjj + AVjjj 



dp 

(M 



E 



'^, cancels out: 

tfW JV W W V 



CJf u 



f 




(52) 



Remembering that, in the approximation under consideration, the vertex functions, O, 
reduce to unity, this total contribution has an attractive character. It just cancels the 
contribution which arises from removing the energy dependence from Eqs. (1, 4, 10) and 
corresponds to the renormalized potential V, and the associated wave function, 0. 
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Figure 11: Selected contributions to the two-body interaction due to two boson exchange 
(crossed-type diagrams). Apart from a common factor l/{uJiUJf(uJi +ujf)), they involve a factor 
l/uJiUf (a,b). 



a) b) 

Figure 12: Example of time-ordered two boson exchange contributions with negative energies 
for the constituents (Z-type diagrams). These contributions, which should be incorporated in a 
full description of the two-body interaction, don't interfer with those considered in this work. 

This result is important in the sense that it indicates that usual single boson exchange 
potentials, that are used in the Schrodinger equation and are described by Yukawa poten- 
tials, already account, in a hidden way, for part of the two boson exchange contribution, 
including crossed diagrams. By the same token, they also account for the renormalization 
and recoil effects related with the energy dependence of the genuine interaction, Ve, but 
this may be true only at the lowest order in ■\. 

Neglecting some non-locality, a physical interpretation of the above contribution to the 
two-body interaction may be obtained by going back to its expression in configuration 
space: 

AVi + AVn + AVni = Vi(r) Vo{r), (53) 

where Vo(r) and Vi{r) have been defined previously (Eqs. 6, 9). It corresponds to the 
interaction between the two constituents, Vo(r), multiplied by the boson in flight proba- 
bility, which is given by Vi{r). It is not a surprise therefore if it cancels the normalization 
correction to the interaction potential, |{V^i(^), Vo{r)}, which was corresponding to miss- 
ing this contribution. Obviously, as already mentioned, the cancellation holds provided 
that the two constituent particles in the intermediate state with a boson in flight have 
the same spin, isospin or angular momentum as the initial or final states. A little bit of 
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caution is nevertheless in order in identifying Eq. (53) with the contribution of Fig. 2, 
which at first sight could be written as the interaction of the two constituent particles, 
Vo(r), multiplied by the probability to have a boson in flight. The argument misses the 
point that, in order to get Vo(r) in Eq. (53), all relative times corresponding to an emis- 
sion of a boson by one constituent and its absorption by the other have to be considered. 
This cannot be achieved by crossed diagrams alone as the time at which the absorption 
takes place (in Fig. 10a for instance) is necessarily bounded by the time at which the first 
emitted boson is absorbed. This is the place where the contribution of the box diagram 
(Fig. 9a) is relevant as it precisely provides the missing time ordered contribution. In 
view of this, one can imagine that the cancellation of the norm correction to the potential 
together with that one involving the interaction of the two constituents while one boson 
is in flight (part of Fig. 2) can be generalized to an undetermined number of bosons 
in flight, within the conditions under which Eq. (53) has been derived (locality of the 
interaction, spin- and isospin-less particles, no Z diagram). A schematic proof is given in 
the appendix. 
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Figure 13: Time ordered diagrams showing the relationship of two boson exchange contributions 
considered here (bottom part) with recoil-norm contributions associated to the interaction with 
an external probe (represented by a thin line in the top part). The different diagrams, a, b and 
c, correspond to crossed, box and renormalization contributions. The double bar on the boson 
line in the diagram c) reminds that only the part contributing to the norm has to be retained. 
The bottom diagrams are obtained by bringing in the top ones the line 3 in coincidence with 
the line 1. 

The above cancellation has been obtained in the case of constituents with equal masses. 
There is no difficulty to check that it also holds with unequal masses within the same 
approximations. In the present approach, this is a first step in allowing one to recover in 
particular the usual Coulomb interaction of an electron (light particle) in the field of a 
nucleus (heavy particle). 

The cancellation of contributions given by Eq. (52) with those due to the "renormal- 
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ization" of both the interaction and the wave function, Eq. (17), reminds the norm- recoil 
cancellation discussed in the literature in the 70 's. This one however concerns the interac- 
tion of a two-body system (in the simplest case) with an external probe and thus involves 
at least three bodies. We examined this situation and indeed find some relationship. The 
difference with the present case is that one of the two bosons exchanged between the two 
bodies in Figs. 9, 10, 11 is attached to a third body, the electron for instance for an 
electromagnetic probe, see Fig. 13. Apart from a possible factor related to the number 
of constituents, there is a one to one correspondence. The demonstration is achieved by 
noticing that all time ordered diagrams are summed up while the two bosons attached 
to the same constituent line interact independently with the two other lines, with the 
consequence that the corresponding vertices necessarily commute. 



4.3 Massless neutral and spinless bosons 

The logarithmically divergent character of the extra contribution to the interaction given 
by Eq. (17) in the limit of a zero mass boson also occurs for the contribution given by 
Eq. (52). Like there, to overcome the difficulty, one has to take into account the binding 
energy in the boson propagators appearing in Eqs. (49-51). A priori, the corrections 
differ from those appearing in Eq. (27), preventing contributions of the two neutral boson 
exchanges to cancel that one due to the renormalization of the wave function. However, 
we saw that the dominant correction arising from this renormalization, of relative order 
a log -, was unaffected by the precise value of the coefficient to be inserted in front of this 
binding energy (2 in Eq. (27)). As a result, this particular correction cancels out when 
considering both the renormalization of the wave function and the two boson exchange 
contribution. 

Remaining corrections to the interaction may therefore be of relative order a. At first 
sight, it is not clear whether these ones should also cancel out from the consideration 
of the above contributions alone (assuming they should cancel at all as expected from 
relativistic calculations where the correction to the dominant order is often believed to 
be of the order a^ rather than a). The problem is that the long range of the force makes 
it difficult to discuss the contributions order by order, the binding energy introduced to 
remove the logarithmic divergence in Eq. (27) accounting for some of them. However, on 
the basis of the result presented in the appendix and assuming we have not missed some 
mathematical sublety, we expect that these corrections of relative order a with respect 
to the dominant term should vanish. It is obviously out of question to discuss here the 
next correction of relative order a^, which would also suppose to precise the framework 
(equation) to be used as well as to account for other radiative corrections. 



4.4 Residual contribution for bosons with spin or charge 

In the case of the exchange of bosons carrying some charge (isospin or color) or some spin, 
or coupling to the spin of constituents, the vertex functions, O, in Eqs. (49-51) do not 
commute. Some residual contribution then remains when adding the contribution of two 
boson exchanges, Eqs. (49-51), and the contribution arising from the renormalization, Eq. 
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(17). It occurs for bosons currently referred to as the pion in the field of nuclear physics 
or the gluon in QCD. Such contributions are generally ignored. Part of them may be 
accounted for in a phenomenological approach, by fitting parameters of a potential to NN 
scattering cross sections for instance. While one cannot underestimate the power of this 
procedure, it should be noticed that the resulting two pion contribution has certainly a 
very complicated structure and it is not sure that it can be easily approximated by single 
meson exchanges. One can thus imagine that these last approaches, by trying to reproduce 
fine contributions which their simplicity does not permit in principle, introduce some bias 
in other sectors of the interaction. With this respect, we notice that the Paris model of 
the NN interaction is the only one where the two pion exchange contribution is correctly 
accounted for [0. It however suffers from some drawbacks, such as neglecting terms of 
the order (p/m,)^, which now turn out to be relevant in determining a parametrization of 
the model p3|. The two pion exchange contribution was also considered in the full Bonn 



model ||T2[ , but its energy dependence and the problems this feature are raising have led to 
its discarding. It is finally worthwhile to mention the description of the Nucleon Nucleon 
system from the Bethe-Salpeter equation [^ or from other relativistic approaches ||37 . 



Based on single meson exchange, they incorporate contributions discussed in the first part 
of this work (Sect. 2), but neglect contributions corresponding to the crossed diagrams 
of Figs. 10 and 11, thus introducing another kind of bias. 

We tentatively examined contributions due to the above residual two meson exchange 
contribution in the case of the NN interaction. The idea is to get an order of magnitude 
and, consequently, to know whether one has to worry about them. In this aim, we neglect 
the contribution of the box diagram, Eq. (49) and a contribution of the same magnitude 
for the crossed diagram. The terms left out are of the order ujiUJf/{uJi + ujf), where 
uji and uof refer to the on-mass shell energies of the two exchanged mesons. This is at 
most equal to 1/4 (cjj = ujf), which gives some estimate of what is neglected. Due to 
possible destructive interferences, actual corrections may be larger however. With the 
above approximation, the expression of the two meson exchange contributions given by 
Eqs. (49-51) becomes identical to that of the norm correction to the energy independent 
interaction, Eq. (17), except for a different ordering of the spin and isospin operators. 
Furthermore, we concentrated on the contribution of isospin 1 mesons, which in particular 
contain the pion whose contribution to the NN interaction is known to be important, and 
pick up the part whose spin structure is the same as for the term 1(^5 K)}- The residual 
contribution under consideration thus reads 

= 4(fI.r^)^{fi,K.}, (54) 

where the first and second terms in the squared bracket respectively arise from the renor- 
malization of the interaction (generalization of Eq. (15)) and the crossed diagrams (gen- 
eralization of Eq. (53)). The 'symbol indicates that only the isospin 1 meson contribution 
to the full interaction is considered, the corresponding product of isospin matrices being 
factored out. 

The relative effect of the two terms is immediately obtained. For the deuteron, of which 
isospin is zero, they are respectively proportional to 9 and 3. This evidences a construc- 
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tive interference, contrary to the neutral spinless meson case where the interference was 
destructive and complete. Amazingly, the second term in the bracket on the r.h.s. of Eq. 
(54) involves the interaction of two nucleons in ^5*0, ^Pq,-- states, which may look surpris- 
ing. One can convince oneself however that nothing is wrong as this state is coupled to a 
meson in flight which, beside some angular momentum, carries isospin 1, allowing one to 
get the total expected isospin, 0. 
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Figure 14: Effect of a residual NN interaction due to isospin 1 mesons: contributions to the 
S and D wave components of the deuteron state shown in the figure are calculated perturbatively. 

The effect of the residual interaction, Eq. (54), can be seen in Fig. 14 for both the S and 
D wave components of the deuteron state (Bonn-QB model). Not much is seen at very low 
momentum for this interaction model. For the S wave, this may be accidental as a large 
effect was obtained with the Bonn-Q model, comparable to the total correction obtained 
previously from the term |{Vi, Vq} (see Figs. 4 and 5). At intermediate momenta, the 
zero at 2 fm~^ is washed out in both cases while at large momenta the contribution blows 
up. It is not clear whether this is a real effect or the result of picking up a part of the 
potential in a domain where what comes from the exchange of mesons with different 
isospin is not well determined. 

Using approximations similar to those made in getting Eq. (54), one can also get an 
expression for the residual interaction due to the gluon exchange (Coulomb part). It reads 



yqq 



(a;.a;)(a;.a;) - a; (ai.a;) x; 



sym. 2 



6{X,.X2)-{V,,Vo}, 



(55) 



where the A matrices act in the color space. The " symbol now indicates that only the 
Coulomb part of the interaction should be considered. Using Eq. (34), the above correc- 
tion amounts to renormalize the Coulomb part of the gluon exchange by a factor roughly 
equal to 0.5 for a^ = 0.5. The size of the effect certainly requires a more careful and com- 
plete examination, especially in view of the fact that it reduces the role of the Coulomb 
part of the one gluon exchange in explaining the meson and baryon spectroscopy. It does 
not support the result of some phenomenological studies that are favoring larger values 
of as, in particular in relation with the position of radial excitations [^, |5^ . 
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Notice that the residual interaction, Eq. (55), involves a log- corrections, thus in- 
dicating that they don't cancel in all cases as the example of the Coulomb interaction 
may a priori suggest. These corrections have some relationship with the well known log 
corrections that account for the renormalization of the weak interaction of quarks due to 
QCD effects. 



5 Conclusion 

In this work, we have studied the relationship of energy dependent and energy independent 
two-body interaction models. The energy dependence of the first ones, based on a field 
theory approach, has obviously some theoretical support but, as is well known, is also 
a source of difficulties. The energy independent models are partly motivated by their 
empirical success, starting with the success of the Coulomb interaction in describing the 
dominant part of the electromagnetic interaction. 

Our study has concentrated on the effect of possible corrections to wave functions. 
Essentially, we showed that effects accounted for by an energy dependent model could be 
incorporated to some extent in an energy independent approach. The resulting interaction 
is different from the usual ones, of the Yukawa or Coulomb type. These ones turn out 
to be renormalized by the probability that the system under consideration be in a two- 
body component, excluding that part involving bosons in flight. Far from raising some 
problem, we, on the contrary, believe that this result provides a useful intermediate step 
in the derivation of effective two-body interaction models to be employed in an energy 
independent scheme. Indeed, beside the above contribution, we expect in any case some 
contribution due to two crossed boson exchanges. For neutral spinless bosons, this one 
cancels the contribution due to the renormalization of the interaction together with an 
extra smaller two boson box contribution. This result, which holds at the zeroth order in 
the inverse of the constituent masses, possibly different and finite, allows one to recover 
the well known Yukawa or Coulomb interactions (at least for some part in the last case). 



It confirms what has been obtained by other methods for these interactions [40, 41 



|3|, Q, perhaps providing a more intuitive explanation. It is not however clear whether 
the conditions required to obtain the cancellation are identical in all cases, authorizing 
some scepticism about the result [^. From present results, from those quoted above [Q 



and from results to come concerning bound states |^T|, one can thus make the following 
conjecture for neutral spinless bosons. The description of low energy states obtained from 
solving the Bethe-Salpeter equation with a single boson exchange interaction as well as 
that obtained with an energy dependent interaction (including on the light front) has 
a strong relationship to a description obtained in an energy independent picture with a 
renormalized Yukawa or Coulomb interaction. Physics described by a Yukawa or Coulomb 
interaction in an energy independent scheme rather has a strong relationship with that 
obtained from the Bethe-Salpeter equation or an energy dependent scheme together with 
an interaction including crossed boson exchange beside single boson exchange. In the case 
of bosons with charge or spin, some residual contribution is left, which by no means seems 
to be negligible for hadronic systems of current interest in nuclear or particle physics. 
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While contributions related to the energy dependence have to be considered and cannot 
be separated from these other ones allowing one to recover the Coulomb interaction for 
instance, what has to be done in practice with them depends on whether one looks at the 
two-body interaction from a theoretical or a phenomenological point of view. In the first 
case, they are part of a series of contributions that have to be determined with more or 
less accuracy depending on the domain and especially on the size of the coupling. In the 
other case, as in the NN strong interaction one, a large body of data has allowed one to 
determine phenomenological interaction models, accounting for a large part of the new 
contributions. Dealing with them requires some care, beyond that one mentioned in the 
literature when employing wave functions issued from an energy dependent model fitted 
to NN scattering data [113], 
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In the energy independent scheme and for the simplest case, this supposes to cancel the 
new contribution with another one equal to it but of opposite sign: nothing is changed with 
respect to the original model. If one considers that the new contribution has a genuine 
character, what we believe, then a change of the parameters determining the shorter range 
part of the models has to be performed. In the energy dependent models, the same can 
be done but in any case a readjustment of the parameters is required. Furthermore one 
has to keep in mind that wave functions in both schemes are not equivalent and cannot be 
used on the same footing as input as well as ouput. Thus, the zeroth order unperturbed 
wave function in an energy dependent scheme should not be taken as that one in an 
energy independent scheme as one may expect at first sight |jl^, |^, [10| , but rather as the 
product of this one by the operator, (1 + Vi)~2, where Vi is related to the derivative of 
the interaction with respect to the energy, Vi = —-^■ 

The present work has been essentially devoted to the first order contributions issued 
from the energy dependence of field-theory based interaction models. In the case of a 
phenomenological approach to the NN strong interaction, they turn out to be largely 
included in these models, as already mentioned. They have an elusive character for some 
part and should not provide much effect in analyzing physical results at low or moderate 
energies. In view of the present work, the large effects found in refs [|1^, thus appear 
as resulting from accounting roughly for twice the same contribution. Beyond first order 
corrections, some effect in relation with boson in flight should show up. They are likely 
to require studies more refined than those presented here. 
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A Summation of all contributions to the interaction 
at the zeroth order in — 

m 

We here discuss the summation of a subset of diagrams contributing to the two body 
interaction in the non-relativistic hmit. Our intent is to show that the cancellation of sec- 
ond order contributions arising mainly from renormalization and crossed boson exchange 
has a somewhat general character for neutral spinless bosons, leaving the single boson 
exchange contribution as the only effective one. Before entering into the details, let's 
mention that this result has been extended to three bosons. Beside all irreducible boson 
exchange diagrams, it supposes to take into account corrections of the order Vi in the 
renormalized interaction, Eq. (13), as well as second order contributions in the expansion 
of the energy denominators given by Eq. (3). 



The developments follow those presented for instance in ref . [E^ , but with a difference as 



to eliminating the elementary interaction of the constituents with the bosons. This one will 
be carried only on the part of the Hamiltonian containing the boson mass term, Ho{b), and 
the interaction term, Hi{c,b). The residual interaction of constituents with the bosons, 
generated by the transformation of the kinetic Hamiltonian of the constituents, Ho{c), 
will be of the first order in the inverse of their mass, implying second order contributions 
for the interaction between the constituents, hence the announced result. 
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The schematic Hamiltonian we start from may be written as: 

H{c, b) = Ho{c) + Ho{b) + Hj{c, b), (56) 

with 

p k pS ^^ 

where the destruction operators, Cp- and 6g, (and the cooresponding creation oprators), 
respectively refer to the constituents and the bosons, whose exchange is responsible for 
the interaction between constituents. 

We now make the substitution 

[c+, c , 6+, b] = e^(^^'^ '^^'^) [C+, C , B+, B] e-^(^^'^ '^■''^), (57) 

where S{C^ , C , B^ , B) is an antihermitic operator given by : 

The above transformation, which allows one to express the bare degrees of freedom (c, 
b) in terms of the effective ones (C, B), is a unitary one, leaving commutation relations 
unchanged. It is no more than a change of basis in the Fock space. With the above 
substitution, the Hamiltonian is written: 

if (c, b) = e^H{C, B) e-^ = H{C, B) + [5, H{C, B)] + i [5, [S, H{C, B)]] + ... (59) 
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Concentrating on that part involving the boson mass and interaction terms in H{c,b), 
one can show that the choice made in Eq. (58) allows one to remove the interaction term 
linear in the coupling g. For this purpose, one uses the relation: 

Hi{C,B) + [S,Ho{B)]=(}. (60) 

Taking into account this relation, the next contribution of the second order in the coupling 
g involves the double commutator, [S, [S,Ho{B)]]. There are two terms arising from the 
second and third terms in Eq. (59), with the factors —1 and | respectively. The resulting 
contribution is nothing but the usual two-body interaction between the constituents: 

H2body{C) = -2 H ^ ^^+k'^P^i-k'^P'- (^^) 



p,]? ,k 



'k 



Beyond the double commutator, it can be checked that the multiple commutators of S 
with Hq{B) vanish, leaving the two-body interaction, Eq. (61), as the only contribution 
from the boson mass and interaction terms of the original Hamiltonian. It is identical to 
the single boson exchange potential in the instantaneous approximation. Notice that to 
obtain this result, S has to commute with H2body{C), which is only fulfilled for spin- and 
charge-less bosons. While the above contribution is quite simple, the contribution arising 
from the mass term of the constituents, Ho{c), is not. At the lowest order in the coupling 
g, it produces an interaction term given by: 

Hj(C,B) = y^E=CtrCp(B+,-BA (^IZ3±R] _ (62) 

As seen from the above expression, the effect the interaction may produce is of the order 

2 

of the difference in the energies of the constituents, namely — , which means an effect 

of the order (— )^ for the interaction of the constituents. This is a typical relativistic 
correction of the order (-)^, expected in any case. As shown in the text of the paper in 
a limited case, the various corrections to the effective two-body interaction beyond the 
instantaneous approximation therefore cancel out at the zeroth order in — for spin- and 
charge-less bosons. 
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